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The Hubbard model in the Heisenberg limit is studied by Kondo-lattice theory. The Kondo 
temperature 7k or /cbTk, which is an energy scale of low-energy local quantum spin fluctuations, 
is enhanced by the resonating valence bond (RVB) mechanism, so that 7k — Tmf/ (2D), 
where Tmf is the Neel temperature in the mean-field approximation of the corresponding 
Heisenberg model and D is the spatial dimensionality. Electrons certainly behave as localized 
spins at T> 7k, but they are still itinerant at T <C 7k unless an antiferromagnetic complete 
gap opens. When the Neel temperature Tk is so high that 7k 3> Tmf/ (2D), magnetism is 
prototypic local-moment magnetism. When 7k is so low that 7k >C Tmf / (2D) because of 
low dimensionality or frustration, magnetism is itinerant-electron magnetism of an almost spin 
liquid, i.e., a normal Fermi liquid or a Tomonaga-Luttinger liquid in which the spectral weight 
of single-particle excitations is almost vanishing. The spin susceptibility has a temperature and 
wave-number dependence characteristic of itinerant-electron magnetism. This type of itinerant- 
electron magnetism must also be possible in the Heisenberg model. 

KEYWORDS: Hubbard model, Heisenberg model, itinerant-electron magnetism, local-moment magnetism, 
Kondo lattice, resonating valence bond, Fermi liquid, Tomonaga-Luttinger liquid, spin liquid, 
Mott insulator 



1. Introduction 

The problem of moment formation or magnetism is 
one of the most important issues in condensed-matter 
physics. 1,2 ) There are two types of magnetism, local- 
moment magnetism and itinerant-electron magnetism. 
According to this nature of magnetism, there are two 
types of effective Hamiltonians, a localized-spin model 
such as the Heisenberg model and an itinerant-electron 
model such as the Hubbard model. It is really a long- 
standing issue how to construct a unified theory that 
can explain the two types of magnetism. When electrons 
are localized, an itinerant-electron model is reduced to a 
localized-spin model. First of all, the mechanism of elec- 
tron localization should be elucidated. 

The theory or scenario for electron localization pro- 
posed by Mott and other people, 3 ) such as the Mott in- 
sulator or the Mott- Hubbard metal-insulator transition, 
is controversial. Consider the Hubbard model with the 
on-site repulsion U and, e.g., only the nearest- neighbor 
transfer integral t. When U/\t\ = +oo and N = L, where 
N is the number of electrons and L is that of unit cells, all 
the unit cells are singly occupied and there is no empty 
or double occupancy, so that no electron can be itiner- 
ant. This ground state is a prototypic Mott insulator. 
It is a highly degenerate state such that it has an elec- 
tron of arbitrary spin at each unit cell or the entropy is 
fcs In 2 per unit cell. The third law of thermodynamics is 
broken in the prototypic Mott insulator. Then, one may 
suspect that, if the ground state is a Mott insulator for a 
finite U, the third law must also be broken in it; and the 
other may suspect that, since the third law must hold, 
the scenario itself is doubtful. 

In the Heisenberg limit, the Hubbard model is reduced 
to the Heisenberg model with the superexchange inter- 
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action J = —4t 2 /U. When J = 0, the ground state is 
certainly the prototypic Mott insulator. Fazekas and An- 
derson 4 ) proposed a stabilization mechanism of a spin 
liquid by the formation of a local but itinerant singlet or 
a resonating valence bond (RVB) on each pair of nearest 
neighbors; this spin liquid is called an RVB spin liquid 
and this stabilization mechanism is called an RVB mech- 
anism. When J < 0, the prototypic Mott insulator is 
unstable against an RVB spin liquid. 

In order to avoid confusion, we tentatively restrict our 
discussion to the ground state within the constrained 
Hilbcrt subspace where no order parameter exists. Ei- 
ther of a Mott insulator and a spin liquid is an insulator. 
They should be distinguished from each other. In this 
paper, only an insulator where the third law is broken is 
called a Mott insulator. If the third law holds in an insu- 
lator, it is Wilson's insulator or a spin liquid according 
to the spectrum of low-energy spin excitations. 

According to Hubbard's theory, 5,6 ) the Hubbard gap 
opens and the band splits into the upper and lower Hub- 
bard bands. When U > 2z\t\ at least, where z is the coor- 
dination number of nearest neighbors, the Hubbard gap 
is a complete gap and the ground state is a Mott insu- 
lator. However, the virtual process allowing empty and 
double occupancies is possible, unless \t\/U = 0. The 
superexchange interaction arises from the virtual pro- 
cess. 7 ) Since the Fock-type term of the superexchange 
interaction stabilizes a local singlet on each pair of near- 
est neighbors, it eventually stabilizes a singlet ground 
state or a normal Fermi liquid 8 ) rather than a Mott in- 
sulator; this stabilization mechanism is simply the RVB 
mechanism. In a previous paper, 9 ) it was proved based 
on Kondo-lattice theory 10-12 ) that the ground state can- 
not be a Mott insulator, unless U/\t\ = +00 and N = L. 
When U is finite, the Hubbard gap is a pseudo-gap. 

Hubbard's theory is within the single-site approxima- 
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tion (SSA). If all the single-site terms are considered in 
a theory, it is simply the supreme SSA (S 3 A) theory. 
The S 3 A is rigorous for infinite dimensions, 13 ~ 16 ) but 
within the constrained Hilbert subspace. According to 
the Kondo-lattice theory, 10 ~ 12 ) the S 3 A is mapped to a 
problem of self-consistently determining and solving the 
Anderson model, 10 ~ 12 ) which is an effective Hamiltonians 
for the Kondo effect. The S 3 A can also be formulated 
as the dynamical mean- field theory (DMFT) 17 ~ 20 ) or the 
dynamical coherent potential approximation (DCPA). 21 ' 
It was also proved in the previous paper 9 -* that, in the 
S 3 A, i.e., even if the RVB mechanism is not considered, 
the ground state is a normal Fermi liquid. The density of 
state has a three-peak structure, a quasi-particle band 
between the upper and lower Hubbard bands, which 
corresponds to the three-peak structure in the Ander- 
son model, the Kondo peak between two sub-peaks. The 
quasi-particle band is simply what is predicted by a com- 
bination of Gutzwiller's theory, 22-24 ) which is also within 
the SSA, and the Fermi-liquid theory. 25,26 ) 

One may suspect that the superexchange interaction 
cannot work in a metal since the formation of localized 
spins is assumed in the well-known derivation of the su- 
perexchange interaction. 7 -* In field theory, the superex- 
change interaction arises from the virtual exchange of a 
pair excitation of an electron and a hole in the lower 
and upper Hubbard bands. 27 ~ 29 ) The superexchange in- 
teraction works in not only a spin liquid but also a metal 
where the Hubbard gap is a pseudo-gap. 

Multi-site terms are perturbatively considered by the 
Kondo-lattice theory. 9 ~ 12 ) A normal Fermi liquid in the 
S 3 A, which is stabilized by the Kondo effect, is further 
stabilized by the RVB mechanism. The stabilized Fermi 
liquid is a normal Fermi liquid, which is totally frus- 
trated. 30 ' The Fermi liquid is a relevant unperturbed state 
to study a normal or anomalous Fermi liquid in the con- 
strained or whole Hilbert space and an eventual true 
ground state in the whole Hilbert space. The Kondo tem- 
perature Ik is an energy scale of low-energy local quan- 
tum spin fluctuations in a normal or anomalous Fermi 
liquid or an ordered state. Since the entropy is about 
k-B In 2 per unit cell at T ^ Tk, a high temperature phase 
at T ^> Tk is simply a Mott insulator. Electrons behaves 
as localized spins at T > Tk- The Kondo-lattice the- 
ory can treat not only itinerant electrons at T <C Tk 
but also localized spins at T > Tk- It is interesting to 
study the Heisenbcrg limit of the Hubbard model by the 
Kondo-lattice theory in order to elucidate the mechanism 
of moment formation and the nature of itinerant-electron 
magnetism and local-moment magnetism. 

One of the purposes of this paper is to show that 
the Kondo-lattice theory based on the Hubbard model 
can explain local-moment magnetism in the Heisenbcrg 
model. Then, it will be shown that, in lower dimen- 
sions, itinerant-electron magnetism is possible even in 
the Heisenbcrg limit. This paper is organized as follows. 
Preliminary is given in §2. The Heisenbcrg limit is stud- 
ied in §3. Discussion is given in §4. Conclusion is given 
in §5. An equality is proved in Appendix. 



2. Preliminary 

2.1 Two effective Hamiltonians 

2.1.1 Hubbard model with an electron reservoir 

In the grand canonical ensemble, the averaged elec- 
tron number is irrational, in general, and is a continu- 
ous function of the chemical potential. An electron reser- 
voir should be explicitly considered to ensure this crucial 
property. In this paper, for simplicity, we assume that 
both the Hubbard model and the reservoir are on a hy- 
percubic lattice in D dimensions. 

The total Hamiltonian is composed of four terms: 

H = H d + H c + V d - C - /i (M d + M c ) . (2.1) 
The first term is the Hubbard model: 

n d = e d J2 ™«x - -j= Y^ d L d J* + U Y1 n ^ n ^ ( 2 - 2 ) 

id v (ij)ry * 

where e d is the band center, d ia and di a are creation and 
annihilation operators of an electron with spin a at the 
ith unit cell, n^ = d ia di ai —t d /\>D is the transfer inte- 
gral between nearest neighbors (ij), and U is the on-site 
repulsion. For convenience, the transfer integral includes 
dimensionality T>, 13 ~ 16 ) so that an effective bandwidth is 
0(| £d|). The second term stands for the reservoir: 

t 



H c 



£ c J2 C l* Ci ° ~ ~7^ E C \a C 3<" ( 2 - 3 ) 

ia V (ij)a 



where every notation is conventional as in Eq. (2.2). The 
third term is an infinitesimally small hybridization be- 
tween the Hubbard model and the reservoir: 



V d - 



X 



Y [ v i d la C i" + v * c L d ic 



(2.4) 



where Vi is a hybridization matrix element at the ith unit 
cell; Vi may be zero or nonzero at a unit cell. For con- 
venience, an infinitesimally small but nonzero numerical 
constant, 



A = 0+, 




(2.5) 


is introduced. In the last term, 






Afd = ^2 d L d ia, K~- 


~ Z_v c ia C iai 


(2.6) 



and n is the chemical potential. The existence of the last 
term means that the Hubbard model and the reservoir 
are embedded in another electron reservoir. 
Throughout of this paper, we assume that 



A* = £d 



1 
2 U 



(2.7) 



Then, (Afd) = (■A/'c) = L, where (• • • ) stands for the 
statistical average and L is the number of unit cells. Only 
the just half-filling case is studied in this paper. In order 
to recover the translational symmetry, for convenience, 
we consider averaged quantities over an ensemble of i>j's, 
which is denoted as {w}, under an assumption such that 
{v} is totally random, i.e., ((«»)) = and ((^w*)) = 

\v \ 2 Sij, where ((•••)) stands for the ensemble average over 
{v}. The periodic boundary condition is assumed. The 
thermodynamic limit of L — > +oo is also assumed. 
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2.1.2 Heisenberg model with a thermal reservoir 

The Heisenberg limit is defined by U/\td\ ~ ► +00 with 
the constraint that 

J = -U 2 d /U, (2.8) 

is kept constant, together with the half-filling condition 
(2.7). In this limit, the Hamiltonian (2.1) is reduced to 

H' = H' d + H c + Va. -nN'c + Eo. (2.9) 

If the Hilbert space is constrained to a subspace where no 
empty or double occupancy exists such that n^+riii = 1, 

S 4 = i]T<x^U T ,, (2.10) 

tt' 

where cr = (o x , o~y, o~ z ) is the Pauli matrix, is simply a 
localized spin at the ith unit cell. The first term is the 
Heisenberg model: 

2D 2^> 

m 



K = -^Y J {s i -s j ). (2.11) 



The second term "H c is given by Eq. (2.3). The third term 
is an exchange interaction between the Heisenberg model 
and the reservoir: 

%-c = -A 2 EE X '( S ' ' ^'R^, (2.12) 

i tt' 

where Ki — —2\vi\ 2 /U. The interaction V' d _ c should be 
relevant in the Heisenberg limit, i.e., 

{K t ) = -2W/U, (2.13) 

should be kept nonzero and finite. This condition should 
be included in the definition of the Heisenberg limit. The 
fourth term is the chemical potential term for the reser- 
voir. The last term is a constant: E = L(—U/2 + J/2). 
In the mean-field approximation, the static susceptibil- 
ity of the Heisenberg model (2.11) obeys the Curie- Weiss 
(CW) law for any wave number q: 

Xs(0,q) = l/[kBT-J B (q)/A], (2.14) 

where 



J s (q) = (2J/VD) l p(q), 



where 



<p(q) 



1 D 

— Y^ cos {q v a) 



(2.15) 



(2.16) 



where a is the lattice constant. Throughout this paper, 
the susceptibility is defined in such a way that the con- 
ventional factor g 2 /ig/4 is not included, where g is the 
g factor and lx-q is the Bohr magneton. According to 
Eq. (2.14), the Neel temperature is T MF = -J/(2fc B ), 
and the ordering wave number of the Neel state is Ql = 
(±7r/a, ±7r/o, • • • , ±7r/a). It will be shown in §3.3.1 that 
this prototype of local-moment magnetism can be repro- 
duced by the Kondo-lattice theory. 

2.2 Kondo-lattice theory 

2.2.1 Kondo effect in the Hubbard model 

When U = 0, the ensemble averaged single-particle 
Green function for electrons in the Hubbard model is 



given, in the wave number representation, by 

1 



«G CT (ie„,fe))) 



(2-17) 



ie„ + fx - e d - E(k) - T d (k r , 

where e n = (2n + l)i:k-QT, with n being an integer, is a 
fcrmionic energy, 

E(k) = -2tMk), (2.18) 

with <p(k) defined by Eq. (2.16), and 

r,(k„) = AV^: 



L ^ ie r . 



\x 



2tMk) 



(2.19) 



In the following part of this paper, the notation ((•••)) 
is not explicitly shown and, e.g., the ensemble averaged 
Green function is simply called the Green function. Since 
A 2 is infinitesimally small, terms of 0(A 4 ) and higher are 
ignored in T d (ie n ). It is obvious that 



Imr d (e + 10) < 0, 



(2.20) 



for |e| < 2\/Z)|£ c |, where T^e + iO) is the analytical con- 
tinuation of T d (ie„). The Fermi surface is a hyper-surface 
composed of such fcp's that satisfy 



f(k F ) = 0. 



(2.21) 



Within the Hilbert subspace where no order parameter 
exists, the Green function for nonzero U is given by 



G da (\e n ,k) = l/[is n + ix-€d- E(k) 

- E CT (ie„,fe) - r d (i£ n )]. 



(2.22) 



where ^^(ien, k) is the self-energy. The self-energy is de- 
composed into the single-site £ CT (ie n ) and the multi-site 
AS (J (i£„,fc): T, a (ie n ,k)_= S .(i£ n ) + AS .(i£ n ,fe). To cal- 
culate the single-site E CT (i£„) is mapped to a problem 
of sclf-consistcntly determining the Anderson model and 
solving it. In general, the Anderson model is character- 
ized by the depth of the level of localized electrons mea- 
sured from the chemical potential, the hybridization en- 
ergy between the localized level and the conduction band, 
and the on-site repulsion. The mapping condition is so 
simple that the depth and the on-site repulsion should 
be equal to e d — fx and U, respectively, of the Hubbard 
model; and the hybridization energy, which is denoted 
by A(e), should be determined to satisfy 



Gcr(i£„) = i?<j(i£„), 



(2.23) 



where 

G a {\.E n ) 



i£V, 



• ix- td- S CT (ie r 



de 



, A(e') 



1£„ - £' 

(2.24) 
is the Green function for the Anderson model and 



Rai}£ n ) = {1/ L)' s y\G da {ie n ,k), 



k 



(2.25) 



is the site-diagonal Green function for the Hubbard 
model. Then, it follows from Eq. (2.23) that 



A(e) = Im S a ( £ + iO) + l/R a {e + iO) 



(2.26) 
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As proved in the previous paper, 9 ' in general 
A(e) > -Imr d (e + iO). 



(2.27) 



According to Eqs. (2.20) and (2.27), A(e) > for |e| < 
2\t c \^/D. Since A(0) > simply means that the Fermi 
surface exists in the conduction band of the Anderson 
model, it is proved that the ground state of the Anderson 
model is a normal Fermi liquid, even if the ground state of 
the Hubbard model is anomalous; the single-site S cr (ie n ) 
is exactly equal to the self-energy for the normal Fermi 
liquid in the Anderson model, even if the multi-site self- 
energy is anomalous. 9 -' 

Since the Anderson model is determined at each T, it 
includes T as a parameter. The parameter T is denoted 
as X" to distinguish it from T itself. In this paper, the 
Kondo temperature 7k , which is an energy scale of low- 
energy local quantum spin fluctuations, is defined by 



[*.(0)] T=OK = l/«teTK(T , ) ) 



(2.28) 



as a function of T", where x s (0) is the susceptibility of the 
Anderson model determined at T', and T — K in the 
subscript means that the Anderson model is embedded in 
a reservoir of T = K. Thus, T K (T') or T K (T) depends 
on T in general. If no confusion arises, 7k (7 1 ) is simply 
denoted as 7k in this paper. Since A 2 > 0, 7k cannot be 
the absolute zero Kelvin, i.e., 7k > K definitely. 

The self-energy for the Anderson model in the presence 
infinitesimally small magnetic field 77 is expanded in such 
a way that 31 ' 32 ) 



f 



1 



S CT (e + iO) = -U+(l- (p-y)e + a(l - 4> s )-gp B H 



, s 2 + (nk B T) 2 

l <Pw 7— 

kb 7k 



(2.29) 



for |e| <C k B Tj<_ and T <C 7k, where 7 , <j> a , and <fi w are 
positive constants. Since charge fluctuations are almost 
completely suppressed in the Heisenberg limit, the Wil- 
son ratio is two: 31-33 - 1 



(2.30) 



W s = s /0 7 = 2. 
According to the mapping condition (2.23), 

p{e) = --ImGJe + iO) = --ImRJe + iO). (2.31) 

Because of the particle- hole symmetry, RcG CT (+iO) = 
Rei? CT (+iO) = 0. Then, it follows that 34 ' 

tt/9(0)A(0) = 1. (2.32) 

If p(0) — > or the ground state of the Hubbard model 
is approaching an insulating or spin-liquid phase, A(0) 
of the mapped Anderson model is diverging according to 
Eq. (2.32); this property is crucial in §3.1 of this paper. 
When the expansion (2.29) is used, the Green function 
(2.22) is written in such a way that 

^ 2 ^ (irk B T) 2 



G { d % + i0,k) = 



Y e - E(k) + i 



7 £ 



k B T K 



A£ ff (e + i0,fc)-r d (e + i0) 



(2.33) 



be approximately used at \s\ < k B T^ and T < 7k- The 
Green function (2.33) is called a coherent part of the 
Green function; it describes only single-particle excita- 
tions in the vicinity of the chemical potential. 

Because of i(j> w e 2 /k B T^, no complete gap can open 
even if AS CT (e + iO, k) diverges for any k as e — > 0. The 
ground state cannot be a Mott insulator nor an insulator 
proposed by Lieb and Wu, 35 -* in which the opening of a 
complete gap is expected. 

Hubbard's theory or the Hubbard approximation is 
only relevant at high energies \e\ ^> /cbTk, 5,6 * 1 i.e., the 
Green function in the Hubbard approximation can de- 
scribe only single-particle excitations in the upper or 
lower Hubbard band. In this paper, such a Green func- 
tion that is only relevant at |e| 3> k B T^ or |e| > k B T^ is 
called an incoherent part of the Green function. 

2.2.2 Intersite exchange interaction 

The irreducible polarization function in spin channels 
ir s (icoi,q) is also decomposed into the single-site 7r s (iw/) 
and the multi-site Air s (iw; , q) such that -k s (iu>i , q) = 
TTsiiui) + A7r s (iwz, q), where u>i — 2lirk B T 1 with I being 
an integer, is a bosonic energy. The susceptibilities of the 
Anderson and Hubbard models are given by 



Xs{m) = 27r 8 (iwj)/[l - Un a (iui)], 



(2.34) 



and 



X s (\ui,q) = 27r 8 (iwj,q)/[l - Uir s (iuJi,q)], (2.35) 

respectively. Since either of x s (iw/) and \si}^i,q) is 
nonzero and finite, Utt s (iuji) — > 1 and Uir s (iu>i , q) — > 
1 should be satisfied in the Heisenberg limit, so that 
ATT s {iuji,q) = 0{l/U 2 ). Then, it follows that 



X s (iuji,q) = 



Xsfal) 



l-(l/4)I e (w,q)x,QuiV 



(2.36) 



where 



I s (iuji,q) = 2U 2 Air s (iuji,q). 



(2.37) 

Here, terms of O [1/Ux s (iuji)] and higher are ignored, 
because they vanish in the Heisenberg limit. 

Equation (2.36) is in full accord with a physical picture 
of a Kondo lattice such that local spin fluctuations at 
different unit cells interact with each other by an intersite 
exchange interaction. The intersite exchange interaction 
is simply 7 s (iw;, q). It is decomposed into three terms: 

I s (iuJi, q) = J s (q) + J Q (icoi,q) - 4A(iwj, g). (2.38) 

The first term J s (q) is the superexchange interac- 
tion, 27 ~ 29 ) which is defined by Eq. (2.15). The second 
term Jq{\uji,q) is an intersite exchange interaction that 
arises from the virtual exchange of a low-energy pair ex- 
citation of an electron and a hole. When the three-point 
single-site reducible and irreducible vertex function in 
spin channels are denoted by A s (ie„ + kj;, ie„; iw;) and 
A s (ft„ + iu>i, ie„; hoi), respectively, it follows that 



A s (ie„ + iuj h ie n ;iuJi) 



X s (ie n + iuJi,ie n ;iui) 
1 - Un s (iu>i) 



(2.39) 



This is accurate at |e| -C ^bTk and T <C 7k and can 
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In the Heisenberg limit, 



2A s (ie„ + ie ni ie n ;iuji) 



UxsKm) 

According to the Ward relation, 36 ) 

d 

dhz 



(2.40) 



A 8 (ie n ,ie n ;0) = l- — &|.(ie„) - S4.(i£„) , (2.41) 



where hz = ff/^B-ff- Thus, it follows from Eqs. (2.29), 
(2.40), and (2.41) that 

A s (i£„,i£„;iu; ( ) = 24> a /[Uxa0ui)], (2.42) 

for \e n \ <C ksTK- When Eq. (2.42) is approximately used 
as the single-site vertex function for \e n \ < ^b^k an d 
l w (| ^5 fcB^kj Jq^uQ) is calculated to be 

J Q (iwi,q) = [4W?/g(wj\[PQwi,q) - FoQwi)] , (2.43) 

where W s — 2 is the Wilson ratio, 

P(m,q) = -^Yj 2 7 G%(ie n +m,k+q)G%(is n ,k), 



and 



L ^^ 

nh<7 



P {m) 



l^P(i^,g). 



(2.44) 



(2.45) 



The single-site term is subtracted in Eq. (2.43) to avoid 
double counting; G^(ie n ,fe) is the coherent part of the 
Green function; and, for convenience, P(iwi, q) is defined 
in such a way that it includes <j> 2 . The third term A(iwj, q) 
is the sum of the whole remaining terms; it corresponds 
to the mode-mode coupling term in the self-consistent 
renormalization (SCR) theory of spin fluctuations. 37 ) 

When Eq. (2.42) is used, an exchange interaction me- 
diated by intersite spin fluctuations is calculated to be 

\ 2 M 0; ™i)U 2 [XaOui, <Z) " Xs(i^)] = 4> 2 s I* s (icj u q), 

(2.46) 

where the single-site term is subtracted to avoid double 
counting and 



IZQuuQ) 



I s {iui,q) 



(2.47) 



l-(l/4)I a (iui,q)x.Qui) 

A perturbative theory in terms of I s (iwj , q) can treat the 
exchange interaction (2.46); I s (iu>i,q) is enhanced into 
I*(iu>i, q) by intersite spin fluctuations and, according to 
Eq. (2.46), <p s is an effective vertex function. The pertur- 
bative theory is simply the Kondo-lattice theory and is 
also a 1/D expansion theory. 

3. Physical Properties of the Heisenberg Limit 

3.1 Fermi surface of an RVB almost spin liquid 

One of the most crucial terms in I s (iwj , q) is the su- 
perexchange interaction J s (q). The first order term in 
J s (q) is the Fock-type term or the RVB term: 

k B T 



A4 RVB )(fc) 



L ^ 

npa' 

xG'^ ) ,(i£„,p)e 



^Ij s (fe-p)((T 



ie„0 + 



(3.1) 



where G d( ^(is n ,p) is defined by Eq. (2.33) and two </ 
arc included as two effective vertex functions. When 



J s (k 



23 

~~D 



y [cos(k u a) cos(p„a) 



+ sin(Ava) sin^a)] , 
is made us of, the RVB term is simply given by 
A£( RVB )(fc) = (3/4)^W 2 E D (23/DMk), 
where 



~u = 07^- Yl v{k)G { ^{ie n , k) 

nk 



L 



(3.2) 
(3.3) 
(3.4) 



For convenience, So is defined in such way that it in- 
cludes </> 7 . When the RVB term is only considered in the 

multi-site self-energy, G d J. (e + io, fc) is simply given by 

G { £ (e + iO, fc) = 1/ [0 7 £ + 2t*Mk) + iO] , (3.5) 

where the life-time term proportional to 4> w and r^E+iO) 
are ignored and 

t d = t d -(3/4)j>. ( W 2 Z D (3/D). 



Then, it follows that 

H D = (1/L)5>(fc)/ [£(*)] 



where 



and 



C(fc) - -2(t* d /^Mk), 



/(£) = l/[e £ /( fcBT ) + l]. 



(3.6) 

(3.7) 

(3.8) 
(3.9) 



The density of state for the coherent part or for |e-| < \t* d \ 
is given by 



P (c) {e) = 



xV E*[ e -«*)]■ 



j i L k 



(3.10) 



An effective bandwidth of p^ c \e) is about \t d \/4>-y. When 
MSItSlMr. 

(a D /\t d \)[l + 0(e 2 )], D^2 



P {c) (e) = 



(a 2 /\t d \)[\n\t* d /(^e)\+0(e°)], D = 2 



(3.11) 

where ap is a positive numerical constant; ap = 0(1). 
The expansion coefficients <f> s and </> 7 , which are single- 
site properties, should be consistently calculated with the 
RVB term. In principle, this problem should be solved 
by self-consistcntly determining and solving the Ander- 
son model. However, we take a different approach to this 
problem. On the basis of Eq. (2.32), we can assume that 



A(£) 



_/ \t d \/(7ra D ) |e| < |*5|/(2^) 
0, |£| > |*5|/(2^) 



(3.12) 



for the Anderson model. Since the probability of empty 
or double occupancy in the Hubbard model, which is 
approximately given by (t d /D)/U 2 provided that the 
virtual processes allowing empty and double occupan- 
cies are explicitly considered, i.e., the RVB mecha- 
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nism is considered, 34, 38 ' should be equal to that in 
the Anderson model, which is approximately given by 
A(O)(|i2|/0 7 )/C/ 2 , an approximate relation such that 

tj/(DU 2 )^A(0)\t* d \/(^U 2 ), (3.13) 

should hold. It follows from Eqs. (3.6), (3.12), and (3.13) 
that 7 ~ a D TTU/[3 2 W*Z 2 D (J/D)]. Then, 

4> S =2^ = 0(DU/\J\), \t* d \=0(U), (3.14) 



and 



p( c '(0) = 1/[ttA(0)] = 0(1/U). (3.15) 



6 7 and |i^| are diverging and 



In the Heisenberg limit, 

p( c >(s) is vanishing, so that the ground state within the 
constrained Hilbcrt subspace is an almost spin liquid. 34 ' 
When D /2, Tk can be evaluated from the Fermi- 
liquid relation for the Anderson model: 31 ' 32 ' 

Xs (0) = 2<^ c '(0), (3.16) 

at T = K. It follows from Eqs. (2.28) and (3.16) that 

P (c) (0) = 1/(2MbT k ). (3.17) 

It follows from Eqs. (3.6), (3.11), and (3.17) that 

k B T K ~ |t;|/(# 7 ) ~ 3S D |J|/(4D). (3.18) 

When T = K, H^'s are calculated to be Si = 1/n = 
0.31831 • • • , E 2 = 2V2/tt 2 = 0.28658 • • • , and E^ = 
l/(2-v/7r) = 0.283095 • • • , i.e., E D ~ 1/3 for any D. When 
D ^ 2, it follows that 



T K ~ T M f/(2£>), 



(3.19) 



where Tmf = |«T|/(2A:b) is the Neel temperature in the 
mean-field approximation of the Heisenberg model. 

When D = 2, x s (0) has a logarithmic singularity due 
to that of p( c ' (e) . The Anderson model is determined at 
each T or T', as discussed above. When the life-time term 
in E CT (e + iO) is considered in the mapping process at T", 
the singularity is suppressed, so that 

[U0)] T=OK * (2a 2 a /|t2|) ln|(ty^)/(fc B T')|. (3.20) 
Then, Tk has a logarithmic T dependence such that 

T K (T) ~ (T M F/4)/ln|(i*/0 7 )/(fc B T)|. (3.21) 

Unless T is extremely low, Eq. (3.19) is relevant even for 
T> = 2. 

Note that Ed given by Eq. (3.4) or (3.7) is a decreasing 
function of T. Since E D -> +0 as T -> +oo K, T K (T) ->• 
+0 K as T ->• Too K. However, when T < T MF /(2T>), 
T < T K (T) is satisfied for any £>, even for D = 2. 

The multi-site self-energy A£ CT (i£ n ,fc) is decomposed 
into the RVB term and the other term 5E CT (ie n , fe): 

AE ff (ie„,fe) = A^ VB \k)+S^(ie n ,k). (3.22) 

The RVB term does not depend on ie n , so that it is quite 
normal. If 5T, a (ie n ,k) is ignored, therefore, the ground 
state is a normal Fermi liquid even in the Heisenberg 
limit. Although the Fermi liquid is an almost spin liq- 
uid, the dispersion relation £(fe) of quasi-particles is well 
defined and the Fermi surface can be unambiguously de- 
fined. The Fermi surface is exactly the same as that for 



U = 0, which is the hyper-surface defined by f(k-p) = 0. 
Since the perturbative scheme in terms of I s (iu>i, q) is 
in parallel to the conventional one in terms of U, it is 
straightforward to calculate 5Y, a (ie n ,k). According to a 
preliminary study, anomalous terms exist in dT, a (ie n ,k) 
for D = 1, which implies that the ground state is a 
Tomonaga-Luttinger (TL) liquid in one dimension. If the 
Fermi surface is defined by (p(kp) — 0, it shows a per- 
fect nesting for any D. Then, electrons may behave as 
a TL liquid in any D dimensions, in particular, in two 
dimensions. 39 ' Even if the ground state is a TL liquid, 
5T, a (e + iO, k) is continuos at e = and its imaginary 
part vanishes at s = 0. Then, 5Y, a (+iO, fep) = for any 
&f that satisfies ip(kp) — because of the particle-hole 
symmetry Within the preliminary study, no anomalous 
term can be found for any D such that the Fermi surface 
cannot be defined because of it. If the Fermi surface can 
be defined, which is quite probable, it must be exactly 
the same as that for U = 0. 

3.2 Metallic conductivity at T = OK 

According to the Kubo formula, 40 ' the conductivity is 
given by 

ctxx(w) - (fi/iw) [K xx {w + iO) - K xx (0)] , (3.23) 

where K xx (lo + iO) is the analytical continuation of 

1 f/3 



K xx (iu>i) 



La 



D 



dTe lulT (e TH j x e- TH j x ), (3.24) 



where ft — l/(k B T) and j x is the x component of the 
current operator defined by 



3 = ^2j(k)n ka , 



ka 



where 



and 



j(k) = -(e/h)(d/dk)[-2t d tp(k)], 



rik 



, = (i/L)5>*-<*-*>4U 



JCT: 



(3.25) 



(3.26) 



(3.27) 



where Ri is the position of the ith unit cell. Equa- 
tion (3.24) is also described as 

K xx (hj t ) = TlTrfe n -(i^), (3-28) 



h 2 Da D ~ 2 



where 



1 f 13 

n x x(i^) = t y^sin(fc x q) sm(p x a) / dre lUlT 

L k P Jo 



^(e^nkae-^np,,,) 



(3.29) 



First, we consider only the RVB term in the multi-site 
self-energy The current vertex j(k) should be consis- 
tently renormalized with the RVB term to satisfy the 
Ward relation. 36 ' The ladder vertex first order in J s (q) 
corresponds to the RVB term, which is first order in 
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J s (q)- When Eq. (3.2) is used, it is easy to show that 
1 2n xx (iwi) 



n X x(iui) 



<f& 1 + (3J/2D)W s 2 7r xx (iu; ( ) 



(3.30) 



where 



7 rri 

■K xx {iwi) = -4>l — — y^ sin 2 {k x a)G^l (ie n , fc) 

nfc 

xG^(ie n + iwi,fe). (3.31) 

For convenience, 7T xx (iu;() is defined in such a way that 
it includes 6 2 . When 



dn. xx (\ui) 
d(iwj) 



d7r £C£C (iw/)/(i(iu;/) 



^7 [l + (3J/2£>)W 2 7T xx (iwi) 
is used, the static conductivity is calculated to be 



(3.32) 



Rea xx (0) 



e 2 (2t d 



S xx (0) 



h? Da D ~ 2 



l + (3J/2D)W*n X x(0) 



(3.33) 



where 
Sxx(0) = 4 lim Re — [tt xx {lj + iO) - tt xx (0)] . (3.34) 

A2, w->0 IUJ 

When Eq. (3.5) is used for G d J.(ie n ,k), it follows that 



Kxx{0) = j X! s ' m2 ( k xa) 



df(e) 



ds 



(3.35) 



e=|(fe) 



where £(fc) is defined by Eq. (3.8), and 

<? xx (0) = ^ V S in 2 (fc x a) f +X> de \-M 
ttL^ J_ oc ds 



x ImG^(£ + iO,fe) 



(3.36) 



When T = K, Eq. (3.36) is reduced to 



21i 



S xx (0) = —-y^s'm 2 (k x a) 



Im- 



1 



■2t*Mk) - r d (+i0). 



(3.37) 

where the reservoir term r<j(+i0) is explicitly shown. 
This is approximately calculated to be 

+oo r -. i2 



S xx (0) 



n\t 



ds 



h 



Im 



1 



e-ilmr d (+i0) 



2\t* d lmT d {+iQ)\ 



(3.38) 



According to Eq. (A-6) in Appendix, it follows that 



td 



— t* 



(3.39) 



l + (3J/2D)W 2 ir xx (0) 
Although t d = 0(VU),t* d = O(U). Then, it follows that 

1*51 



Re<r xx (0) 



2hDa D - 2 Imr d (+i0) 



(3.40) 



Since t* d = 0(U) and Imr d (+i0) = 0(X 2 U), Rea xx (0) = 
O (t/°/A 2 ). The conductivity diverges as A 2 —J- 0. Even 



if i5S -(ie„, fc) is considered in addition to the RVB term, 
this conclusion does not change, as studied below. 

When the renormalized current vertex is denoted by 
./(fc; \s n , ie n + iu>i), it follows that 

K xx (m) = -^Vj.(fc)G^(i£„,fc)G^(i en +i^,fc) 

nkcr 

x J x (fc;ie n ,ie n +iwc). (3-41) 

According to Eq. (3.32), each 7r xx (iwi) independently 
contributes to the iw; linear term of Ii X x(^i)- Accord- 
ing to a similar argument, as regards the iu>i linear term, 
it is sufficient to consider 
k B T 



K 'xx(iui) 



La 1 



J2G^(ie n ,k)Gi c \ie n +iui,k) 



nka 



x J 2 (fc;ie„,ie„), (3.42) 

instead of K xx (}0Ji). According to the Ward relation, 36 ) 

(3.43) 



e d 



J(fc;ie„,ie„) = -t^t E d {k) + AS CT (ie„,fc) 

The relation (3.39) is consistent with Eq. (3.43). Al- 
though |j(fc)| = O(VU), |J(fc;ie„,ie„)| £ ^ - 0(U). 
If the current vertex is consistently renormalized with 
the multi-site self-energy, Rea xx (0) = 0(U°/X 2 ). The 
conductivity at T = K is nonzero and is diverging as 
A 2 — > even in the Heisenberg limit. 41 ) 

3. 3 Two types of magnetism 

3.3.1 Local-moment magnetism, in higher dimensions 

According to Eqs. (2.36) and (2.38), the inverse of the 
static susceptibility is given by 

(3.44) 
Then, we consider a function 7n(<z) defined by 



[l/ X s(0,q)] T=1 



= 0. 



(3.45) 



JT=T N (q) 

The Neel temperature In is simply the highest Tj^(q) as 
a function of q; q that gives the highest 7n is an ordering 
wave number of the Neel state. 

When Wilson's result for the s-d model 33 ) is extended 
to the s-d limit of the Anderson model, it follows that 



1 



X.(0) 



k B (T + 9 K ) 



T<^T K 
T>T K 



(3.46) 



where 9k > and 9k = 0(7k); the definition of 7k in 
this paper is different by a numerical factor from Wilson's 
definition. When T <C 7k, electrons locally behave as a 
normal Fermi liquid. When T 3> 7k, electrons at each 
unit cell behave as a free localized spin. 

First, we consider infinite dimensions (D — > +oo). 
Since T K ^+0KasD^ +oo, l/x s (0) = k B T. Since 
T/Tk — > +oo for any nonzero T, Jq(0, q) is vanishing for 
any q. 42 ' Since D — > +oo, A(0, q) is vanishing for any q. 
Although J s (q) is 0(J/>/~D) and vanishing for almost all 
q's, it is 0(J/D°) and nonzero for particular q's; e.g., 



- J s (0) = J S (Q L ) = -23 = 4fc B 7i 



MF. 



(3.47) 
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where Ql = (±7r/a, ±7r/a, 
easy to show that 



±7r/a). Then, it is quite 



Tk — T] 



MF, 



(3.48) 



the ordering wave number of the Neel state is Ql, 
and the susceptibility at T > Tmf is exactly given by 
Eq. (2.14). In infinite dimensions, the susceptibility of 
the Heiscnbcrg limit exactly obeys the CW law for pro- 
totypic local-moment magnetism. Thus, magnetism in in- 
finite dimensions is prototypic local-moment magnetism. 
In finite dimensions, Tk — Tmf/ (2D). The reduction 
of Tn from Tmf is a higher order effect in 1/D. When 
D >• 1, therefore, it is obvious Tmf 3> Tk and Tn S> Tk. 
When T > T N , i.e., T > T K , Jq(0, q) is small for any 
q. When D » 1, A(0,q) is also small. At T > T N , the 
inverse of the susceptibility is approximately given by 

« H ^ + ^-i J ^ (3 ' 49) 

The Neel temperature is approximately given by 



Tn — Tmf — 9i> 



(3.50) 



where 9 K = O [T MF /(2D)}. When T N > T K , magnetism 
is local-moment magnetism even in finite dimensions. 

3.3.2 Itinerant- electron magnetism in lower dimensions 
Because of critical fluctuations, no order appears at 
T > K in one and two dimensions. 43 ' Then, Tn can 
be low in actual quasi-one and quasi-two dimensional 
magnets; Tn can also be low in actual highly frustrated 
magnets. Here, we assume that Tn <C Tk- When Tn < 
T < T K , Eq. (3.44) is reduced to 



fc B T K -ij s (g)-ij Q (0,<7) 



Xs(0,q) 



A(0,q). 

(3.51) 

This can be approximately used at T < Tk- The T de- 
pendence of l/xs(0, q) arises from Jq(0, q) and A(0, q). 
When 5H a (ie n ,q) is ignored, the density of states for 
quasi-particles is defined by 



(£) = ^E (5 [ £ -^ fc )] = ^p (c) (£), 



(3.52) 



which is well defined even in the Hcisenberg limit. Since 
4W^7x s (0) 2 = 4 2 (fc B T K ) 2 , it follows that 



l - J Q (0, q) ~ 4(fc B T K ) 2 [P(0, g) - P o (0)] . 



The subtraction term is given by 

P (0) = 2 fds 1 fe 2 p*(e 1 ) P *(e 2 ) 



f( ei)-f(s2) _ 
ei-e\ 

The T dependence of -Po(O) is small. 

The homogeneous (q = 0) and staggered (q 
components of P(0, q) are given by 

df(e) 



P(0,0) = 2 dep*(e) 



ds 



and 



P(0,Qz 



cfep* (e) - tanh (fie) 



(3.53) 

(3.54) 

= Ql) 
(3.55) 

(3.56) 



When D y^ 2, p*(e) is almost constant around e ~ 0, as 
shown in Eq. (3.11), it follows that 



P(0,0) ~2p*(0)[l + 0(T 2 )}, 



(3.57) 



and 



P(0, Q L ) * 2p*(0) [ln(T K /T) + 0(T°)] , (3.58) 

When D = 2, p*(s) is logarithmically diverging as s — *■ 0, 
as shown in Eq. (3.11), it follows that 

2a 2 r, 



P(0,0) 



and 



1*31 



In |^/(fc B T) | + O(T ) 



(3.59) 



P(0, Q L ) ~ -^ { [In |0 7 ^/(fc B T)|] 2 + O(T )} . (3.60) 

The T dependence of P(0,Ql) or Jq(0,Ql) is larger 
than that of P(0,0) or Jq(0,0) because of the perfect 
nesting of the Fermi surface. The strength of Jq(0, q) is 
proportional to the bandwidth of £(fc) or fc B TK. 42 ' 

If either of SJ^ a (ie n , k) and A(0, q) can be ignored, the 
staggered susceptibility shows a logarithmic T depen- 
dence due to the perfect nesting of the Fermi surface, i.e., 
it approximately obeys the CW law; this is a mechanism 
of the CW law in antiferromagnetic metals. 44, 45 > Since 
the density of state p*(e) has a sharp peak at e — in 
two dimensions, the homogeneous susceptibility approx- 
imately obeys the CW law; this is a mechanism of the 
CW law in ferromagnetic metals. 46 - ) 

In general, the q dependence of A(0, q) is small. Even 
if (*>£ CT (£ + iO, fe) is considered, the Fermi surface shows a 
perfect nesting for Ql, so that the temperature depen- 
dence of Jq(0, Ql) must be much stronger than that of 
Jq(0, 0). Thus, the qualitative features of Xs(0, Ql) and 
Xs(0,0) when 5Y, a (e + iO, fe) is considered must be the 
same as those when 5E CT (e + iO, k) is ignored; the T de- 
pendence of Xs(0,Ql) is stronger than that of Xs(0, 0). 
Such q dependence is characteristic of itinerant-electron 
magnetism, and it can be distinguished from the q depen- 
dence of local-moment magnetism, where only the Weiss 
constant depends on q. 

3.3.3 Magnetization in the Neel state 

When we follow previous papers, 45 ' 47 ' it is straight- 
forward to derive Landau's free energy below Tn- In two 
dimensions and higher, provided that critical fluctua- 
tions are ignored, staggered magnetization at T < Tn 
and T ~ Tn is given by 



KjJ« Vi-(7YTn) s 



(3.61) 



either when Tn *C Tk or when Tn ^> Tk- 

When Tn 3> Tk, an antiferromagnetic (AF) gap much 
larger than A: b Tk opens at T < Tn- The spectrum of 
single-particle excitations is not well defined at any T. If 
an AF gap is smaller than /s b Tk, which is only possible 
when Tn ^C Tk, the spectrum of single-particle excita- 
tions is still well defined at T <C Tk, even at T < Tn, in 
the Hcisenberg limit. 
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4. Discussion 

The Kondo temperature Tk or fceTk is the energy scale 
of local quantum spin fluctuations. When T 3> Tk, local 
thermal spin fluctuations overcome the quantum ones so 
that electron behave as localized spins. When Tn S> Tk, 
therefore, magnetism is local-moment magnetism. When 
T <C Tk, on the other hand, electrons are itinerant even 
in the Heisenberg limit, in general. An antiferromagnctic 
gap opens at T < Tn- Since the gap is a complete gap 
in the model of this paper, the conductivity vanishes at 
T = K, i.e., electrons are localized at T = K. If 
the gap is much smaller than /crTk, however, electrons 
can still behave as itinerant ones so that magnetism is 
itinerant-electron magnetism. If the gap is much larger 
than fceTK, which is presumably only possible when 
Tn 3> Tk, electrons cannot behave as itinerant ones so 
that magnetism is local-moment magnetism. Whether 
magnetism is itinerant-electron or local-moment mag- 
netism is determined by whether Tn <C Tk or Tn S> Tk- 

It is an interesting issue whether or not the above dis- 
cussion on the Heisenberg limit can be extended to the 
Heisenberg model. The relationship between the Heisen- 
berg and Hubbard models is similar to that between the 
s-d and Anderson models. According to Yosida's the- 
ory, 48 ) the ground state of the s-d model is a singlet. 
Anderson's scaling theory supports the singlet ground 
state. 49 ) Wilson proved that the ground state is a singlet; 
a crossover occurs between a localized spin at T> Tk 
and a local spin liquid at T <C Tk- 33 ) Nozieres proposed a 
Fermi-liquid description of the local spin liquid based on 
the phase shift of conduction electrons due to scattering 
by the local spin liquid. 50 ) Either of the ground states 
in the s-d and Anderson models is a singlet or a nor- 
mal Fermi liquid. Then, essentially the same Fermi-liquid 
analysis as Nozieres's was made based on the Anderson 
model by Yamada and Yosida. 31 ' 32 ) In the Heisenberg 
limit of the Hubbard model, the ground state within the 
constrained Hilbcrt subspace is also a singlet, a normal 
Fermi liquid or a Tomonaga-Luttinger (TL) liquid. The 
Green function for the reservoir is given by 

G ca (e + iO, k) = 1/ [e + 2t cV (k) - T c (e + iO)] . (4.1) 

When |e| < fc B T K , it follows that 



r c (e + iO) 
A 2 






i 



-e(fc)-<5S CT (e + iO,fe)/0 7 

(4.2) 



Since \v\ 2 /4> 7 = 0{U°), T c (e + iO)/A 2 is well defined 
even in the Heisenberg limit. The quantity r c (e + iO) de- 
pends on vanishing fcrmionic excitations in the Hubbard 
model, and it corresponds to the phase shift discussed by 
Nozieres. 50 ^ A vanishing fcrmionic spectrum exists in the 
s-d limit of the Anderson model or the Heisenberg limit 
of the Hubbard model. It is reasonable that there is a 
trace or vestige of the vanishing fcrmionic spectrum even 
in the conduction band of the s-d model or the reservoir 
for the Heisenberg model. 

The charge susceptibility is exactly zero in the s-d 
model but is nonzero in the Anderson model. The con- 
ductivity is exactly zero in the Heisenberg model but is 



nonzero, finite or infinite, in the Hubbard model. These 
differences are obvious because the local electron num- 
ber at each unit cell is a conserved quantity in the s-d 
and Heisenberg models but is not in the Anderson and 
Hubbard models, i.e., local gauge symmetry exists in the 
s-d and Heisenberg models but it does not in the An- 
derson and Hubbard models. On the other hand, local 
gauge symmetry can never be spontaneously broken or 
recovered. 51 ' The recovery of local gauge symmetry in the 
s-d or Heisenberg model is caused by a non- spontaneous 
process, i.e. by constraining the Hilbert space within the 
subspace where no empty or double occupancy is allowed 
at a localized spin site or each unit call. Because of this 
peculiar nature of local gauge symmetry, it is never a rel- 
evant symmetry to classify or distinguish a phase from 
the other. Thus, the adiabatic continuation 52 ) holds be- 
tween a spin liquid in the s-d model and an electron 
liquid in the Anderson model, i.e., low-energy physical 
properties of a spin liquid in the s-d model can be de- 
scribed as those of a normal Fermi liquid in the Anderson 
model. According to the scaling theory by Abrahams et 
al., 41 ) there is no mobility edge between itinerant and lo- 
calized states in a disordered system or there is no min- 
imum metallic conductivity, which means that there is 
no qualitative difference between a metal and an insula- 
tor. Thus, the adiabatic continuation can hold between a 
metal and an insulator, in general; e.g., it holds between 
Wilson's insulator and a doped metal. The adiabatic 
continuation must also hold between an electron liquid 
in the Hubbard model and a spin liquid in the Heisen- 
berg model such that low-energy physical properties of a 
spin liquid in the Heisenberg model can be described as 
those of an electron liquid in the Hubbard model. Thus, 
it is reasonable to speculate that itinerant- electron type 
of magnetism is also possible in the Heisenberg model. 
Since Tk — Tmf/(2T>), where Tmf is the Neel tempera- 
ture in the mean-field approximation for the correspond- 
ing Heisenberg model and D is the spatial dimensional- 
ity, whether magnetism is itinerant- electron type of mag- 
netism or local-moment magnetism must be determined 
by whether T N < T MF /(2D) or T N > T MF /(2D). 

It has been proposed that, in one dimension, low- 
energy spin excitations in the Heisenberg and XXY mod- 
els can be mapped to those of a TL liquid. 53 ~ 56 ) This pro- 
posal is simply a proposal that magnetism in the Heisen- 
berg and XXY models must be itinerant-electron type of 
magnetism. According to the Kondo-lattice theory, the 
spectrum of low-energy pair excitations in the Heisen- 
berg limit arc determined by the imaginary part of the 
exchange interaction Jq(uj + iO,q) or P(w + iO,q). When 
the RVB term is considered but (JE^e + iO, q) is ignored, 
the energy a; pa j r (g) of a pair excitation lies in 

63i|J| sm(qa) < w pair (g) < 12Si| J\ sm(qa/2), (4.3) 

as a function of wave number q. The lower limit corre- 
sponds to des Cloizeaux-Pearson's mode in the Heisen- 
berg model, 57 ) whose dispersion relation is given by 



Udcp(q) = (2| J\/tt) sin(ga). 



(4.4) 



When 5Y, a (e + iO, q) is ignored, Si = 1/ir, as shown in 
§3.1. If Si = 1/ir, the lower limit is three times as large 
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as cj(jcp(9)- The imaginary part of S'E a (s + iO, q) reduces 
Ed. The energy dependence of (5£ CT (e + iO, q) also reduces 
S^. It is interesting to study how w pa i r (q) is corrected or 
modified when 8Y, a (e ± iO, q) is considered. 

In the Heiscnbcrg model on the square lattice, the Neel 
state is only stabilized at T = K: Tn = +0 K. The stag- 
gered susceptibility Xs(0,Ql) is diverging as T — >• K. 
On the other hand, the homogeneous Xs(0, 0) shows a 
peak around T MF ^ | J|/fc B - 58 ~ 60) This feature, x s (0,g) 
having qualitatively different T dependences for differ- 
ent q's, is characteristic of itinerant-electron magnetism. 
For example, a similar result was also obtained by a 
conventional perturbative theory in terms of U based 
on an itinerant electron model. 61 ) When we follow pre- 
vious papers, 45 ' 46 ' 61 ^ it is easy to show that, since the 
Fermi surface shows a perfect nesting and the density of 
state p*(s) has a sharp peak at the chemical potential, 
A(0, q) has a large T dependence such that it suppresses 
the Curie- Weiss T dependence for any q; the q depen- 
dence of A(0, q) is small. The reduction of Tn and the 
peak around Tmf — \J\/kB hi Xs(0, 0) can also be ex- 
plained in terms of the mode-mode coupling A(0, q) by 
the Kondo-lattice theory based on the Hubbard model. 
Furthermore, magnetization at T = K is as small as 
0.3 /zb per unit cell. 62 -* This small number implies that, 
in the corresponding Heisenberg limit of the Hubbard 
model, an antiferromagnetic gap is smaller than fceTk so 
that magnetism is itinerant-electron magnetism. Thus, 
we propose that, even in the Heisenberg model in two 
dimensions, magnetism at T <C Tmf/4 should be char- 
acterized as itinerant-electron type of magnetism. 

5. Conclusion 

The Heisenberg limit of the Hubbard model is stud- 
ied by the Kondo-lattice theory. The Kondo temperature 
Tk or fceTK can be still defined as an energy scale of 
low-energy local quantum spin fluctuations even in the 
Heisenberg limit. It is enhanced by the resonating va- 
lence bond (RVB) mechanism, so that Tk ~ Tmf/(2-D), 
where Tmf is the Neel temperature in the mean-field ap- 
proximation of the corresponding Heisenberg model and 
D is the spatial dimensionality. When T 3> Tk, elec- 
trons behaves as localized spins and the entropy is about 
fee In 2 per unit cell. Thus, a high temperature phase at 
T 3> Tk is a Mott insulator. Within the constrained 
Hilbert subspace where no order parameter exists, how- 
ever, the ground state is not a Mott insulator nor a Lieb- 
Wu insulator either but is an almost spin liquid, i.e., a 
normal Fermi liquid or a Tomonaga-Luttinger liquid in 
which the band of low-energy single-particle excitations 
is almost vanishing; however, the width of the vanishing 
band is 0(A;bTk), the Fermi surface is well defined in 
the vanishing band, and the conductivity at T = K is 
diverging in the vanishing limit of disorder. 

When D 3> 1, the Neel temperature Tn is so high 
that Tn — Tmf but the Kondo temperature is so low 
that Tk <C Tmf- Thus, Tn S> Tk, so that magnetism 
for D ^$> 1 is local-moment magnetism. In particular, 
magnetism in infinite dimensions, where Tk = +0 K 
and Tn = Tmf, is prototypic local-moment magnetism. 
When D is small enough, Tk is so high or Tn is so low 



that Tn <C Tk- Electrons are itinerant at T <C Tk, 
unless an antiferromagnetic complete gap opens in the 
vanishing band. Magnetic properties at Tn < T <C Tk 
and \u>\ <C ^bTk are those of a normal Fermi liquid 
or a Tomonaga-Luttinger liquid, i.e., an itinerant elec- 
tron liquid; e.g., the spin susceptibility has a temperature 
and wave-number dependence characteristic of itinerant- 
electron magnetism. When Tn <C Tk, therefore, mag- 
netism is itinerant-electron magnetism. 

Since local gauge symmetry is never relevant to clas- 
sify a phase from the other, there is no essential differ- 
ence between an electron liquid in the Hubbard model, 
where local gauge symmetry does not exist, and a spin 
liquid in the Heisenberg model, where local gauge sym- 
metry exists, except for charge- channel properties such 
as the conductivity. Thus, itinerant- electron type of mag- 
netism must also be possible even in the Heisenberg 
model. Whether magnetism is itinerant-electron type of 
magnetism or local-moment magnetism must be deter- 
mined by whether T N < T MF /(2T>) or T N > T M f/(2T>). 
Since no Tn exists in one dimension and Tn = +0 K 
in two dimensions, in particular, magnetic properties of 
the Heisenberg model at T ^C Tvif and \u>\ <C ^b^mf in 
one dimension and two dimensions must be able to be 
described as those of an itinerant-electron liquid, i.e., a 
Tomonaga-Luttinger liquid or a normal Fermi liquid. 
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Appendix: Proof of Eq. (3.39) 

Consider a function defined by 

y^>(fc)/[£(fc)-M*] 



a(m* 



L ^ 



cos 



(k ia )f[ak) - u? 



(AT) 



where £(fe) = — 2{t* d /4> 1 )ip{k). Equation (A-l) is also 
given, in the integration form, by 

\f~Dn D r +7r / a f+^/a 

W^T^ir/ rffc i"7 dfc D cos(fcio)/ [£(*)-/**]. 

^ ) J—7r/a J—Tt/a 



(27 



(A-2) 



By partial integration of Eq. (A-2) with respect to fci, it 
follows that 



X(n*) = 2(t* d /4> J )Y(n*), 



(A-3) 



where 



K(/0 



D r+ir/a r+ir/a 

dki ■ ■ ■ dkfj sin (k\a) 

— Tr/a J—-n/a 



(2ir) D 



df(e) 



de 



(A-4) 



e=£(fe)-M* 

This is also given, in the sum form, by 

dm 

L ^r-' de 



y ^*) = 7E sin2 ( fcia ) 



e=£(fe)-ji* 



J. Phys. Soc. Jpn. 



Full Paper 



Fusayoshi J. Ohkawa 11 



1 1 ° 

k v=l 



df(s) 



de 



e=f(fe)-A»* 

(A-5) 



Since Sr; = X(0) and w xx (0) = Y(0), it follows from 
Eq. (A-3) that E D = 2(t* d /4> 7 )Tr xx {0) or 



2 /. 3 7 ~ 2 _ J 



Then, it is straightforward to prove Eq. (3.39). 
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